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ABSTRACT 

Aims. To investigate in detail the morphology at low radio frequencies of the supernova remnant 

(SNR) IC 443 and to accurately establish its radio continuum spectral properties. 
O . 

5_l Methods. We used the VLA in multiple configurations to produce high resolution radio images 

of IC 443 at 74 and 330 MHz. From these data we produced the first sensitive, spatially resolved, 
spectral analysis of the radio emission at long wavelengths. The changes with position in the ra- 
dio spectral index were correlated with data in near infrared (NIR) from 2MASS, in gamma-rays 

^ ■ from VERITAS, and with the molecular 12 CO (7=1-0) line emission. 

IT) • 

Results. The new image at 74 MHz has HPBW = 35", rms=30 mjy beanr 1 and at 330 MHz 
CN ! HPBW = 17" and rms=1.7 mjy beam The integrated flux densities for the whole SNR are 

"^74MHz = 470 ± 51 Jy and Sj^m = 248 ±15 Jy. Improved estimates of the integrated spectrum 
were derived taking into account a turnover to fit the lowest frequency measurements in the lit- 
erature. Combining our measurements with existing data, we derive an integrated spectral index 
a \oMm Hz = ~ 039 ± 001 with a free-free continuum optical depth at 330 MHz t 330 ~ 7 x 10~ 4 
(Tio = 1.07); all measurements above ~ 10 MHz are equally consistent with a power law spec- 
trum. For the pulsar wind nebula associated with the compact source CXOU J061705. 3+222127, 
we calculated S ™^ Hz = 0.23±0.05 Jy, Shmmhz = 0.20±0.04 Jy, andtt^° M ^ z ~ 0.0. Substantial 
variations are observed in spectral index between 74 and 330 MHz across IC 443. The flattest 
spectral components (-0.25 < a < -0.05) coincide with the brightest parts of the SNR along 
the eastern border, with an impressive agreement with ionic lines as observed in the 2MASS J 
and H bands. The diffuse interior of IC 443 has a spectrum steeper than found anywhere in the 

SNR (-0.85 < a < -0.6), while the southern ridge again has a flatter spectrum (a 0.4). At 

the available statistics the VERITAS y-ray emission strikingly matches the CO distribution, but no 
clear evidence is found for a morphological correlation between the TeV distribution and radio 
emission. 

Conclusions. The excellent correspondence between the eastern radio flattest spectrum region 
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and NIR ionic lines strongly suggests that the passage of a fast dissociating J-type shock across 
the interacting molecular cloud dissociated the molecules and ionized the gas. We therefore con- 
clude that thermal absorption at 74 MHz (7-74 up to ~ 0.3) is responsible for the localized spectral 
index flattening observed along the eastern border of IC 443. Towards the interior of IC 443 
the spectrum is consistent with those expected from linear diffusive shock acceleration, while 
the flatter spectrum in the southern ridge is a consequence of the strong shock/molecular cloud 
interaction. 

Key words. ISM:individual objects: IC 443, PWN CXOU J061705.3+222127-ISM:supernova 
remnants-radio continuum:ISM-Infrared:ISM-Gamma rays: ISM-ISM: clouds 



1. Introduction 



The im portance of radio continuum studies of SNRs for und erstanding shock acceleration pro 



cesses ([Reynolds & Ellison 



1992: 



Anderson & Rudnick 



19961) as well as intrinsic and extrinsic 
interactions with ionized gas and the interstellar medium (ISM) ( Dulk & Sleefl975 



Kassim 1989) 



has been long appreciated. However the magnitude of the observable effects are often subtle, re- 
quiring a large leverage arm in frequency space to discern. The lack, until recently, of sufficient 
angular resolution and sensitivity at the lowest frequencies has made progress in such studies slow. 
The advent of high resolution, low frequency observations with instruments like the VLA and 
GMRT are now changing this picture. In particular, VLA sub-arcminute resolution imaging be- 
low 100 MHz has been important for discerning spatially resolved i ntrinsic and extrinsi c thermal 



of sight towards W49B dLacev et al 



1995), thermal 



absorption. Examples now inclu de unshocked ejecta in side Cas A (IKassim et al. 
filaments within the Crab nebula (Bie enhol z et al .1 1997k and ionized gas in the ISM along the line 



Brogan et al.l (12005b spatially resolved the 



12001b . Later on, 

ionized boundary marking the SNR/molecular cloud (MC) interface in 3C 391, suggesting such 
signatures co uld be both common a nd important for delineating elusive SNR/MC interactions. 
Most recently, 



Castelletti et al 



(120071) detected strong 74 MHz free-free absorption at the interface 
between the SNR W44 and the photo dissociation region of a neighbouring HII region. 

In this paper we extend such high resolution, low frequency radio studies to the classic 
SNR IC 443, one of the clearest a priori known cases of a remnant interacting with its cloudy 
surroundings. We present the first, high resolution radio studies of this object, and interpret our 
results in context of recent IR and high energy observations. IC 443 is a particularly attractive SNR 
for low frequency studies both because of its large angular size and because its outer Galaxy lo- 
cation leave it relatively well isolated from the confusing effects normally contaminating studies 
towards inner Galactic SNR/MC complexes. 



1.1. An overview of SNR IC 443 

Observed in the radio domain IC 443 (G189. 1+3.0) is among the larger SNRs in angular size cat- 
alogued in our Galaxy. It consists of two connected, roughly spherical, shells of radio synchrotron 
emission, which are centered at different locations (the names "shell A" and "she ll B" are some 
times used to refer their locations towards the east and west halves of the remnant. 



Braun & Strom 



1986bl) . On its eastern side, IC 443 has a rim-brightened morphology, while towards the western 
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half th e surface brightness is dimmer. A third incomplete and faint shell (called by 



Braun & Strom 



1986bl "shell C") is also evident extending beyond the northeast periphery of the remnant when 



the radio images are displayed with very high contrast. Based on both its morphology and soft X- 
ray spec trum, "shell C" was proposed to be a different SNR called G189.6+3.3, that overlaps with 



IC 443 (lAsaoka & Aschenbach 



19941) . At present it is still questionable if it is a SNR at all, and 



if it is, then it is unclear whether G 189.6+3.3 is interacting with IC 443 (lAsaoka & Aschenbach 
Leahvll2004 



1994 



Lee et alJ 2008) 



At optical wavelengths the appearance of IC 443 resembles that observed in the radio band. 
Spectroscopic and photometric observations of this remnant reveal a complex filamentary structure 
varying in brightness and shape, which is mainly composed of Balmer-li ne emission, strong [Call] , 



[SII], [Nil], [OIII], as well as faint [Fell] and [Felll] emission regions ( Fesen & Kirshner 



1980) 



The emission-line ratios are consistent with those of a well-evolved SNR (Fesen & Kirshne r 1980). 



Kinemat ical consideration s based on optical systemic velocities place the SNR between 0.7 
and 1 .5 kpc dLozinskaval 1 98 11) . although under the assumption of a physical relationship of IC 443 
with t he nearby HI I region Sharpless 249, it might be located slightly more distant between 1 .5 and 
2 kpc ( Fesen|[l984 ). A distance of 1.5 kpc is adopted in most of the works relat ed to IC 443 . 



Dif ferent X-ray obser vations c arried out with the Einstein Observat ory (iPetre et al 



Ginga dWang et al 



1 992 ), ROSAT dAsaoka & Aschenbach 



1994 . ASCA dKawasaki et al 



1988) 



2002) 



and XMM-Newton dTroia et al.1 120061) show a bulk of thermal centrally peaked X-ray emission 



within the radio rim. On the basis of the radio/X-ray morphology and the X-ray properties, several 
author s proposed that IC 4 43 belongs to the class of "mixed-morphology" or "thermal-composite" 
SNRs jRho&Petrelll998l) . 

The compact X-ray source CXOU J061705. 3+222127 is thought to be the remnant of the ex- 
plosion that gave rise to IC 443. This compact relic, discovered using data at high radi o frequencies 



and in X-rays, is embedded within the remnant and is located close to its southern edge dOlbert et al 



20011). Albeit no pulsations were detected from the neutron star, the measure of its transverse ve- 



locity provi des an estimate fo r the age of both the SNR and its associated neutron star, of about 
3 x 10 4 yr (Olbert et al.ll200ll) . An alternative estimate of ~3000 yr was previously proposed by 



Petre et al 



(1988) based on a model of its X-ray emission. 
IC 443 is a prototypical case of a SNR impacting dense interstellar molecular gas. The giant 
molecular cloud, as mapped by CO and HCO + observations, forms a ring in the foreground and 



appea rs to be interacting with the remnant at several locations from north to south (IDickman et al 



1992). This phenomenon, which seems to be responsible for the multiwavelength observational 



picture of the remnant, was unambiguously co nfirmed by the presence of several O H (1720 MHz) 



masers as well as IR cooling lines from H2 dRho et al 



2001 



Hewitt et al 



2006, and references 



therein). Additionally, HI observations show well-defined filamentary struc tures in the north eastern 



region of the remnant, where the brightest optical filaments were located ( Lee et al 



features are interpreted as the recombined neutral gas in an atomic shock dLee et al 



2008). These 



200 



Another interesting characteristic of IC 443 is the presence of associated high e nergy sources. 



The E GRET GeV source (3EG J0617+2238) was detected overlapping the remnant (lEsposito et al 



1996). In addition, TeV gamma-rays from the central region of the remnant, coincident with the 



maser emission, we re recently reported from observations c arried out with M AGIC, VERITAS, and 



FERMI telescopes dAlbert & edl2007 



Acciari et al.ll2009 



Abdo et al.ll2.Q10h . It has been suggested 
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Table 1. Observational summary 



Program 


Observing 


VLA 


Bandwidth 


Integration 


Synthesized 


Largest detectable 




dates 


Config. 


(MHz) 


time (hr) 


beam (arcsec) 


structure (arcsec) 


74 MHz Parameters 


AG0697 


2006 Apr 8, 9, 10 


A 


6.3 


9 


25x20 


800 


AG0735 


2006 Oct 27 


C 


6.3 


8.35 


228x214 


7500 


AG0697 


2007 Oct 21 


AB 


6.3 


6.85 


28x30 


1300 


330 MHz Parameters 


AG0697 


2005 Sep 16 


C 


6.3 


4.7 


62x60 


1800 


AG0697 


2006 Apr 8, 9, 10 


A 


1.5 


9 


6x6 


170 


AG0735 


2006 Oct 27 


C 


1.5 


8.35 


57x53 


1800 


AG0735 


2007 Apr 15, 16 


D 


6.3 


5.2 


238x205 


4200 


AG0697 


2007 Oct 21 


AB 


1.5 


6.85 


22x17 


500 



that the correlation of gamma rays with molecular gas arises from the pion decay of ha dronic cos- 



mic rays generated by t he interaction of the SNR shock with dense molecular material (IHumenskv 



2008 



Abdo et al. 2010). 



2. Observations and data reduction 

The Very Large Array (VLA)Q was used in multiple configurations to observe the large SNR IC 443 
at 74 and 330 MHz. A summary of the observations carried out on various dates between 2005 and 
2007 is given in Table [T] 

In order to help with the radio frequency interference (RFI) excision and mitigate the ef- 
fects of bandwidth smearing all the data were acquired in multi-channel continuum mode (64 
and 16 channels per polarization at 74 and 330 MHz, respectively). At both, 74 and 330 MHz, 
we first performed a bandpass calibration on either 3C 405 (Cygnus A), 3C 147 (0542+498) or 
3C 274 (Virgo A) using publicly available models^. The initial antenna-based gain and phase cor- 
rections at 74 MHz were estimated from observations of either 3C 405 or 3C 274. For the 330 MHz 
data, observations of 3C 147 were sufficient for both, flux density and i nitial phase calibra tion at 



Baars et al 



(1977) scale. 



this frequency. In all cases the absolute flux scale was set according to the j 
The NRAO Astronomical Image Processing Software (AIPS) package was used to process all the 
observations. 

Large fields of view are a general characteristic when observing at long radio wavelengths 
and it is necessary to avoid distortions introduced in the image caused by the non-coplanarity of 
the baselines of the array. To overcome this, we employed a pseudo-three-dimensional Fourier 
inversion as implemented in the AIPS task IMAGR, in which the primary beam area of the 74 and 
330 MHz data, 1 1°.5 and ~ 2°. 5, respectively, is divided into multiple partially overlapping facets 
( Cornwell & Perlevlll992 ). Furthermore, to compensate for time variable ionospheric phase effects 



1 The Very Large Array of the National Radio Astronomy Observatory is a facility of the National Science 
Foundation operated under cooperative agreement by Associated Universities, Inc. 



2 See http://lwa.nrl.navy.mil/tutorial/ to obtain source models in FITS format. 



4 



Castelletti et al.: High-resolution radio study of SNR IC 443 at low radio frequencies 



we performed successive rounds of self-calibration and imaging to the data from each configuration 
separately at 74 and 330 MHz. Angle-invariant self-calibration (as implemented in the AIPS task 
CALIB) is generally inadequate for compensating f or ionospheric effec ts across a large field view, 
particularly for the 74 MHz B and A configuration ( Kassim et al.ll2007 ). An exception is for fields 
dominated by a bright source at the field center, as in the case of IC 443. The final calibrated 
visibility data from all of the observations listed in Table [TJfor each frequency were then combined 
into a single uv data set, after which a further amplitude and phase self-calibration was performed. 
For the concatenated 330 MHz data we employed a multi-scale CLEAN algorithm in AIPS, with 
three different scales sizes. This process is efficient to make high resolution images that are also 
sensitive to extended structures. All the resulting facet images were stitched together into one large 
field using AIPS task FLATN to create a single final image with a synthesized beam of 17 ".27 x 
15". 81, PA=-24°, whic h represents an o r der of m agnitude im provement in angula r resolution over 
the 330 MHz image of 



Braun & Strom (1986b) and that of 



Hewitt et al 



(120061) . The sensitivity 
achieved in our image after correcting for primary beam attenuation is 1 .7 mJy beam -1 . 

Concerning the observations at 74 MHz, we noted that the extended radio emission in the 
field was more properly imaged using t he capability inco rporated in the task IMAGR to switch 



back and forth between the SDI Clean (iSteer et al 



1984) and the usual CLEAN deconvolution 



strategy, depending on the contrast between the brightest residual pixel and the bulk of residual 
pixels after each major cycle. Following this procedure, the final resolution of the first image of 
IC 443 obtained at 74 MHz after combining the data sets described in Tableland including primary 
beam corrections is 36". 50 x 31". 76, PA=34°.71. The rms noise level in the 74 MHz image is 
30 mJy beam -1 . 

Ionospheric refraction and self-calibration normally introduce arbitrary position shifts on low 
frequency images. These shifts are readily corrected by registering against background small di- 
ameter sources whose positions are known from higher frequency observations. We corrected these 
shifts by measuring several small-diameter sources with respect to their known positions from the 
NRAO VLA S ky Survey (NVSS) (the latter has an astrometric accuracy better than 1" in both 



RA. and dec. 



Condon et al 



1998). We determined and corrected for a mean positional differ- 



ence of S .25 + S .06 in R.A. and 5". 15 + 1".00 in dec. at 74 MHz and S .08 ± S .04 in RA. and 
3".71 + 0".42 in dec. at 330 MHz. 



3. Results 

3.1. Low frequency VLA images of SNR IC 443 

In Fig.Q]we present a close-up view of the radio continuum emission at 74 MHz from IC 443. These 
observations provide the first subarcminute image of this remnant created at meter wavelengths. 
The data at 74 MHz are sensitive to smooth structures up to ~ 125' in size, larger than the full 
extent of the SNR, and thus are sensitive to the largest scale structure present in the SNR. 

Figure [2] displays with a resolution of ~ 17" the new VLA image of the ~ 2.5 square degrees 
field of view centered on the remnant observed at 330 MHz^. An enlargement showing the detailed 
330 MHz total intensity morphology of IC 443 is shown in Fig. [3] The combination of data from the 

3 The FITS files of the radio images of IC 443 at 74 and 330 MHz are available online at 
http://www.iafe.uba.ar/snr/SNR.html 
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different VLA array configurations ensures that all scales of the radio emission are well represented 
in the 330 MHz image. 




18m00.0s 6h17m00.0s 16m00.0s 
(J2000) Right Ascension 



Fig. 1. Radio continuum image of IC 443 at 74 MHz. This map has been corrected for the atten- 
uation of the primary beam. The angular resolution of this image was smoothed to a beamsize of 
50". The grayscale is linear ranging from 25 to 250 mJy beam -1 and the contour levels are 72, 
120, 160, 200, 240, and 280 mJy beam -1 . The plus symbol marks the position of the compact 
source CXOU J061705. 3+222127. These observations provide the first subarcminute image of the 
remnant created at meter wavelengths. 



The distinctive morphological characteristics of IC 443, that comprise both the rim-brightened 
eastern side and the breakout into a region of significantly lower density to the west, were 
known from past radio observations. For example, a lower resolution WSRT 327 MHz image 



(Bra un & Stromlll986bl Fig. 5) provides excellent surface brightness sensitivity to the most ex- 
tended structures. Our 330 MHz image (Fig. [3} is quite complementary, with the higher angular 
resolution resolving for the first time at such low frequencies small scale structures covering almost 
the entire extent of the eastern radio shell of the remnant, a region where previous observations 
have indicated t he presence of ionic shocks as well as abundant shocked HI and optical filaments 
( Lee et aljboosl ] t is remarkable that towards the northeast and southeast ends of the eastern shell 
the brightest emission distribution becomes significantly narrower in appearance. An average di- 
ameter of ~ 35', or ~ 15 pc at the assumed distance of 1.5 kpc, is measured in our map for the 
bright radio shell in IC 443 (that, as mentioned in Sect. 1.1 it is referred to in literature as shell 
A, and is depicted in the inset in Fig.|2|. The low surface brightness emission gradually decreases 
outwards to the west forming a more uniform radio shell of about ~ 52' in diameter or ~ 23 pc (the 
so called shell B, Fig. O, as measured on the new 330 MHz view of the remnant. 

Additionally, the 330 MHz image of IC 443 serves to show that even at low radio frequencies 
part of the east outer border of the remnant is quite structured. It is particularly remarkable the 
indented morphology of the bright rim at R.A.= 06 h 18 m 15 s , dec.= +22° 37', with tenuous radio 
synchrotron emission detected ahead of the main shock. Such faint emission is seen as a weak 
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22m00 Os SH20m00.ps "IBmOOOs 16m00 Os 14m00 Os 

(J2000) Right Ascension 



Fig. 2. The image shows the entire field of view of the VLA around the SNR IC 443 at 330 MHz. 
The image displayed includes primary beam correction. The synthesized beam is 17" x 16" with 
a position angle of -24°, and the sensitivity level is 1.7 mJy beam 1 . This is the first high fidelity 
and high resolution view of the emission from IC 443 ever obtained at low radio frequencies. The 
inset with the image of IC 443 at 330 MHz is included to help in the location of the bright and 
weak radio shells referred to in the literature as shells A and B, respectively. 



and irregular radio halo that in our 330 MHz image is notable at a mean level of ~ 4cr extending 
up to 5' ahead of the bright sharp boundary. It is interesting to note that similar weak, diffuse 
emission upstream of the main shock is also observed in the cases of Puppis A and W44, a long 

and 



the sides where the expanding blast wave encountered molecular gas (ICastelletti et al 



2006 



Castelletti et al. 2007). A series of small protrusions are also evident along most of the extension of 
the eastern rim. The largest of these features located at R.A.= 06 h 18 m 20 s , dec.= +22° 38' emerges 
about 2' radially from the border of the remnant. The se structures can be also seen in the image 
of IC 443 at 1420 MHz presented bv lLee et al.l ([2008), who refer to them as "spurs". The authors 
interpret that both, the halo and the "spurs", are likely originated from the interaction of the SNR 
with the surrounding inhomogeneous ambient medium, although an alternative explanation based 
on a physical association between IC 443 and the SNR G189.6+3.3 is not entirely ruled out in 
their work. The low surface brightness emission from G189.6+3.3, positionally coincident with 
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40m00.0s 



S 



22d30m 00.0s 



20m00.0s 



10m00.0s 




18m00.0s 6h17m00.0s ISmOO.Os 

(J2000) Right Ascension 



15m00.0s 



Fig. 3. A color image of the radio continuum emission from IC 443 at 330 MHz constructed us- 
ing multiple-configuration VLA observations. The brightness range covered by the scale is be- 
tween 3 and 35 mJy beam -1 . The final beam size shown at bottom left is 17" x 16" at a posi- 
tion angle of -24°. The noise level is 1.7 mJy beam -1 after primary beam correction. The color 
scale runs between 3 and 45 mJy beam -1 . The white plus sign marks the location of the source 
CXOU J061705. 3+222127. 



IC 443, is clearly detected in our 330 MHz image extending outside the IC 443's eastern shelj^ 
near R.A.= 06 h 18 ra 15 s , dec.= 22° 46'. 

Towards the center of the SNR, in the southern half, the most remarkable feature is a bright 
annular filament placed near R.A.= 06 h 17 m 10 s , dec.= +22° 24', which is immersed in the faint 
and diffuse emission that dominates the inner part of IC 443. Such annular structure forms the 
central part of a more extended feature commonly referred in the literature as the southern sinuous 
ridge. The large southern ridge is defined below dec.~ +22° 35' at 330 MHz and appears to be 
an extension of the southeast border of the remnant. Shocked CO gas with b road lines has been 



detected near this region of the SNR (ICornett et al . 1977 



Dickmanetal. 1992) 



Integrated flux densities estimates for IC 443 were made using the new observations, yielding 
S74 mhz = 470 + 51 Jy and S330 mhz = 248 ± 15 Jy. The quoted values have been corrected for the 
primary beam response and for the contribution of unrelated point sources overlapping the remnant. 



4 The radio continuum emission from G189.6+3.3 is better re presented in all its extension in the low 
resolution image at 327 MHz presented by 
to enhance it. 



Braun & Strom { 1986b), displayed with the appropriate contrast 
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The uncertainties in the measurements account for the statistical errors as well as the selection of 
integration boundaries. 

3.2. The pulsar wind nebula 

The new sensitive image at 330 MHz offers, for the first time, a view of the low frequency coun- 
terpart of the low-luminosity plerionic nebula powered by the source CXOU J061705. 3+222127 
observe d both at higher radio frequencies with the VLA and in X-rays wi th Chandra and XMM- 



Newton dQlbert et al .1 12001 



Bocchino & Bvkov 2001 



v: 



Gaensler et al. 2006). 



A close up image of the pulsar wind nebula (PWN) at 330 MHz is displayed in Fig. [4] At 
74 MHz the combined effect of the lower sensitivity, poorer angular resolution, and flatter in- 
trinsic spectrum, conspire against its visibility at this low frequency. As in X-rays, at low radio 
frequencies the source CXOU J06 1705. 3 +222 127 (marked in Fig. |4]with a plus sign) is placed at 
the apex of the nebula, although a radio counterpart for the point X-ray source is not detected at 
330 MHz. The maximum embedded in the radio nebular emission, with an intensity that peaks up 
to 25 mJy beam -1 , lies about 0'.25 northeastward from CXOU J06 1705. 3 +222 127. The cometary 
shaped nebula extends along its major axis approximately 1 '.25 behind CXOU J06 1705. 3 +222 127 
and 0'.21 ahead of it. The integrated flux density over the entire nebula obtained from our image 
at 330 MHz is S = 0.23 ± 0.05 Jy. In addition we estimat ed the flux densi ty of the PWN 



Lee et al 



(120081) performed by 



330 MHz 

at 1420 MHz on the basis of the image of IC 443 presented by 
combining VLA observations and single dish data taken from the Arecibo Telescope, obtaining 
SrSU = 0-20 + 0.04 Jy. 



22h22m15s 



22m00s 



a 21m45s- 




06h17m10s 17m09s 17m08s 17m07s 17m06s 17m05s 17m04s 
(J2000) Right Ascencion 



Fig. 4. A close-up image of the nebular emission around CXOU J061705. 3+222127 at 330 MHz. 
The position of CXOU J06 1705. 3 +222 127 is shown by the plus sign. The grayscale varies from 19 
to 25 mJy beam -1 . The contour levels on the image are traced at 19, 21, 23, and 25 mJy beam -1 . 



3.3. Radio and optical emission of IC 443 

Figure illustrates, with a spatial resolution better than 20", the morphological comparison be- 
tween the optical and radio emission from IC 443 obtained from the combination in a false color 
image of our new VLA observations at 330 MHz (in red) with data from the Second Palomar 
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Observatory Sky Survey (in green). In this figure, features where both spectral bands overlap are 
shown in yellow. 

The optical emission tracing the low density atomic gas reproduces the east-west asymmetry 
characteristic of the radio total intensity emission as well. The coincidence of the synchrotron- 
enhanced radio emission with very strong optical filaments observed in H a and [SII] lines towards 
the east, suggests that the optical features delineate the position of cooling post-shock ISM gas. 
Furthermore, near infrared emitting gas observed in this portion of the SNR (see below Fig. [9} 
together with abundant neutral gas, correlated both in space and velocity with the optical filaments, 
indicate the presence of radiative shocks propagating into gas of different densities. 

The quality of the new 330 MHz image allows us to identify the close radio/optical correspon- 
dence in most of the small-scale radio structures observed as extensions from the brig ht eastern 
portio n of the shell. This is in good agreement with the behavior previously noticed by 
(2008) using radio continuum observations at 1420 MHz. 



Lee et al 




Fig. 5. A high resolution comparison between radio and optical emission from IC 443 SNR. The 
green corresponds to optical emission from the Second Palomar Observatory Sky Survey, while in 
red the 330 MHz radio emission is shown. The yellow regions are areas where emission in both 
spectral bands overlap. 

The southern radio ridge is also mimicked by a quite faint optical counterpart. Absorption due 
to the molecular gas mainly located in the foreground of the SNR is probably responsible for the 
observed weakness in the optical emission in this part of IC 443. 
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In the breakout region toward the western side of IC 443 only few local radio enhancements, 
immersed in faint diffuse radio emission, are observed at the locations of the optical filaments. 

4. Radio spectral properties of SNR IC 443 

4.1. The spatially resolved spectral index of IC 443 

In this section we analyze the spectral properties of this SNR using the new VLA images presented 
here at 74 and 330 MHz. Estimates of the spectral index variations with position across the remnant 
require imaging the interferometric data of IC 443 at both frequencies using the same uv coverage. 
To perform this, we reconstructed the interferometric images by applying appropriate tapering 
functions to the visibility data at 74 and 330 MHz. To show more clearly the main spectral features, 
avoiding any masking effect from small scale variations, we have chosen a final synthesized beam 
of 70". In addition, to avoid any positional offsets, the images were aligned and interpolated to 
identical projections before calculating spectral indices. 

The dependency of the spectral index between 74 and 330 MHz with the position within IC 443 
was determined through the construction of a spatially resolved spectral index map, which is shown 
in Fig. [6] To produce the spectral map the matched images of IC 443 at both frequencies were 
masked at the 3<x level of their respective noise levels. The error in the determination of the spectral 
index from the map is less than 0.04 in the high flux regions (the east limb of the SNR and some 
interior filaments) and about 0. 1 in the diffuse central emission. The uncertainties increase in the 
weakest emission regions towards the westernmost part of the remnant due to the lower sensitivity 
of the image at 74 MHz in this region (see below). 

Figure [6] shows good agreement between total intensity features and the spatial distribution of 
the spectral index. The 74/330 MHz spectral map reveals for the first time the morphology of a 
very flat spectral component running along the eastern side of IC 443 within which the spectral 
index varies between ~ -0.05 and ~ -0.25. This flattening towards the brightest filaments 
(see also Fig. Q] and Fig. [3} indicates that some thermal absorption is present, and would become 
stronger at lower frequencies. It is remarkable that such spectral behavior has also counterparts in 
the J and H bands as observed by 2MASS (further details for the correlation between local spectral 
index variation across IC 443 and the IR emission are described below in Sect. 5). 

The large-scale diffuse emission in the SNR's interior has a spectrum that is markedly different 
from the eastern part of the remnant with steeper components ranging from a^ = -0.6 to a quite 
steep value of a^P = -0.85, as would be expect ed under the linear diffusive shock acceleration 



1993). The southern sin- 



model for weak shocks with low Mach numbers ([Anderson & Rudnickl l ' 
uous ridge, where most of the interaction with the molecular cloud is taking place, is seen in our 
spectral index map as a region with varying between ~ -0.25 and ~ -0.5. These local spectral 
indices are in good correlation with total intensity features: the brighter synchrotron areas are sys- 
tematically flatter than the other parts of the ridge. In this case the interpretation of the flattening in 
the spectrum is different from the case noted above in the east rim. Here, the spectral behavior might 
be a signature of Fermi shock acceleration at the sites where stronger post compression shock den- 
sities, accompanied by higher local Mach numbers, and/or higher magnetic f i eld strength result due 



to the impact of the SNR blast wave on denser ambient medium (IB ell 
1993). 



1978 



Anderson & Rudnick 
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Fig. 6. Spectral index map constructed using VLA observations at 74 and 330 MHz (matched to 
a common angular resolution of 70")- The 5.5, 15, 26, and 40 mJy beam -1 contours from the 
330 MHz image are included to facilitate the comparison between spectral continuum and total 
power features. To create this map the 74 and 330 MHz images were masked at 3cr. 



Towards the western side of IC 443, which is dominated by diffuse and faint radio emission, the 
spectral indices are mostly steep (o^° - -0.60). Although some flatter features are also observed in 
this region, they are mainly located bordering clipped areas in the map and should not be considered 
as real features owing to the decrease in the signal-to-noise ratio observed in this region in the 
VLA 74 MHz image in comparison with the image at 330 MHz. 

The spectrum calculated for the point-like sources near R.A.= 06 h 17 m 31 s , dec. -+22 ° 45' 48" 



and 06 h 17 m 36 s , +22° 24' 57" confirms their extragalactic nature (IBraun & Strom 



1986a). A sim- 



ilar result is obtained for the point radio source located at the position R.A.~ 06 h 15™ 34.7 s , 
dec- +22° 41' 46". 



4.2. The integrated spectrum of IC 443 

To update the mean spectral index determination for IC 443, we have included the new integrated 
flux densities at 74 and 330 MHz in the extensive list of measurements presented in the literature. 
In Table [2] we list the integrated flux density estimates for the SNR between 10 and 10700 MHz. 
We applied a correction factor over the w ide spectral r ange b etween 408 and 10700 MHz in order 



to place each value on the flux scale of iBaars et al.1 d 1977b . except in some cases for which no 



information was available on the flux value considered for the primary calibrators. 

A plot of the integrated radio continuum spectrum for the SNR IC 443 is shown in Fig. [7] 
Our new integrated flux measurements at 74 and 330 MHz are indicated by filled circle symbols. 
From the spectrum it is evident that the flux densities measured at the lowest radio frequencies, 
and particularly the flux value at 10 MHz, lie below the general trend of the data suggesting the 
presence of thermal absorption along the line of sight. In order to fix the integrated spectral index 
of IC 443 we first use a single power law spectrum to fit the spectrum well down to our lowest 
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Table 2. Integrated flux densities on the SNR IC443 



Frequency Scaled flux References 



(MHz) 


density (Jy) 


10 


400 + 100 <a) 


20 


600 ± 144 (a) 


22 


615 + 75 {a) 


22.25 


535 ± 65 (a) 


22.3 


529 ± 36 (a) 


25 


630 ± 132 <a » 


26.3 


600 ± 48 (a) 


26.7 


561 + 33 (a) 


33.5 


582 ± 37 (a) 


34.5 


440 ± 88 (a) 


38 


650 ± 95 (a) 


38 


730 ± 100 (a) 


38 


460 ± 46 (a) 


38.6 


547 ± 40 (a) 


74 


470±51 (b) 


81.5 


420 ± 63 (a) 


81.5 


470 + 70 (a) 


83 


470 ± 80 (a) 


102 


480 ± 80 (a) 


111 


440 ± 80 (a) 


151 


280 ± 35 (a) 


159 


270 + 40 (a) 


178 


210±42 (a) 


195 


290 ± 45 (a) 


330 


248 ± I5 <b) 


400 


230 ± 34 (a) 


400 


210±31 (a) 


400 


251 ±8 (a) 


408 


289 + 28 



Frequency 
(MHz) 

430 

513 

610 

635 

740 

750 

960 

960 

1000 

1390 

1400 

1400 

1410 

1419 

1420 

1420 

2650 

2700 

3000 

3125 

4170 

5000 

5000 

6640 

8000 

10700 . . . 



Scaled flux 
density (Jy) 



References 
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215 ± 32 ( 
179 ± 18 

164 ± 15 ( 
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177 + 15 
170 ± 20 <c) 
146+ 18 
131 ± 13 
130 ± 13 
160+ 16 (c) 
138 ± 15 
86 ±9 
104 ± 15 
100 ± 15 (c) 
100 ± 15 (c) 
100 ± 15 (c) 
79 ± 11 

85 ± 13 (c) 
70 ± 15 (c) 
90 ± 18 
60 ± 5 (c) 
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(a) No correction to 



Baars etal 



(119771) scale was applied. 
<b) Flux density scale from VLA Calibrator Manual, http://www.aoc.nrao.edu/~gtaylor/csource.html 
(c) The correction factor was not available. 



measurement at 74 MHz, excluding the flux density value at 10 MHz (represented as a solid line in 
Fig. 17). A weighted fit produces a spectral index a = -0.36+0.02 (S v oc •/"). This result agrees very 



Erickson & Mahonev 



(1985 



well within the error limits with the previous estimates presented by 
and references therein). If we consider the lowest frequencies observations, these measurements 
can be fit with a power law plus an exponential turnover using the Eq. Q] (indicated by the dotted 
line in Fig. [7} 
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Sv = 5330 (sso^mJ exp[ - T33 ° (330^) 2 " ] (1) 

Here, 330 MHz is a reference frequency at which an integrated flux density S 330 and an optical 
depth T330 are measured, a represents the non-thermal integrated spectrum, which is assumed to 
be constant throughout the radio band. We have made a weighted fit of the distribution of data 
points over four decades in frequency and find for the whole SNR a single radio spectral index 
a = -0.39 ± 0.01, and an average optical depth 7*330 = (7 ± 1) X 10 4 . The free-free continuum 
optical depth at 10 MHz derived from the relation T10 = T330 [10/330]~ 21 is T10 = 1.07, while the 
optical depth at 74 MHz is T74 = 0.02. Our results indicate that although the absorption becomes 
significantly stronger at 10 MHz this effect is negligible at 74 MHz. In addition, we note that the 
spectral index produced by considering free-free absorption is consistent with that derived with 
a power law. Such a result is not surprising given the results of Sect. 4.1, that clearly indicated 
the presence of thermal absorption, though not at a level to impact the integrated flux at 74 MHz. 
Future measurements at intermediate frequencies (e.g. ~ 30 MHz) are needed to understand if the 
turnover inferred from the 10 MHz measurement is related to the subtle absorption revealed in the 
spectral index analysis presented in Sect. 4. 1 . 
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Fig. 7. Radio continuum spectrum for SNR IC 443 obtained from the flux density values listed in 
Table|2] The filled circle symbols correspond to the new flux density measurements calculated using 
the VLA data at 74 and 330 MHz presented in this work, the open symbols are for ra dio observa 
tions p reviously published and, where possible, brought onto the flux density scale of 



Baarsetal 



d 19771) . Solid line represents the linear fit to the flux density values excluding that at 10 MHz, which 
produces a spectral index a = -0.36 + 0.02 (S oc v°). Dotted line shows a fit to all of the plotted 
values if absorption were present (eq.[T), which yields a spectral index a = -0.39 + 0.01. 



Additionally we re-calculate the spectrum of the pulsar wind nebula in IC 443. To analyze 
the global spectrum of the PWN in IC 443 we combin e the new flux densities estimated at 330 
and 1420 MHz with data taken from loibert et al.l J2OO1I) at 1460, 4860, and 8460 MHz. The radio 
spectrum including all these measurements is shown in Fig. [8] The larger error bars for the lower 
frequency radio data reflect the difficulty in separating the nebular emission from its surroundings. 
A weighted fit to all the flux densities produces a spectral index for the PWN of a = -0.04 + 
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0.05, which is zero within the uncertainty in the fit itself. This result is similar to that obtained by 



Olbert et al 



(2001). 
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Fig. 8. Radio continuum spectrum for the pulsar wind nebula (PWN) around 
CXOU J061705.3+222127 in the SNR IC 443. A weighted fit to the data with a single 
power law function (S oc v°) yields a spectral index a = -0.04 + 0.05. The filled circles symbols 
are from the flux densi ty estimates at 330 and 1420 MHz presented in the current work, while the 
open symbols are from lOlbert et al.l (120011) . 



5. Radio spectral index and the near infrared emission 



Based on our accurate, spatially resolved radio continuum 74/330 MHz spectral map, we have 
investigated the correlation between radio spectral features and the near infrared emission (NIR). 
Figure [9^ shows two 74/330 MHz spectral index contours (traced at a^f = -0.05 and = 
-0.25) enclosing the east rim where we found regions with very flat spectrum in IC 443 superposed 
onto the NIR detected in the J (1.25 nm), H (1.65 //m), an d K s (2.17 fim) bands as taken from the 
Two Micron All Sky Survey (2MASS) jRho et alJboOlb . The 2MASS image of IC 443 shows 
the dramatic contrast in near infrared color between the east rim and the southern portion of the 
remnant. In the color representation of the infrared emission, blue traces the /-band flux, while the 
infrared data in the H and K s are shown in green and red, respectively; white thus enlightens regions 
where all the three IR bands overlap. To facilitate the comparison, Fig. |9p displays the 330 MHz 
continuum image of IC 443 with a grayscale selected to emphasize the brightest radio emission. 

An impressive agreement is observed in location, size and shape between the NIR emission 
detected in the H and J bands and the flattest spectral feature as traced by the a™ contours along 
the eastern edge of IC 443. This correspondence begins in the northernmost part of the remnant and 
extends down to positions near dec.~ +22° 25'. From Fig.|9p it is also notable that the brightest 
ra dio synchrotron e mission perfectly matches the bright emission in the NIR bands. As noticed 
by 



Rho et al 



(120011) . at this site the predominant constituent of the emission in the J and H bands 
is the [Fell] line, with a minor co ntribution from other multi-ionized species like [Nell], [Nelll], 
[Sill], [SIII], etc. lRho et al.l ( 1200 II) proposed a model in which the infrared emission from the ion- 
ized species in the east bright radio limb of IC 443 comes from shattered dust produced by a fast 
dissociating J-type shock. The present accurate comparison between radio spectral indices and IR 
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emission confirms this model. In effect, the passage of a dissociative shock through a molecular 
cloud not only dissociates molecules but also ionize the atoms. Such collisional ionization is re- 
sponsible for the thermal absorbing electrons that produce the peculiar very flat spectrum areas 
observed all along the eastern border of IC 443. This interpretation is in agreement with studies 
based on CO and X-ray o bservations which conclude that the l arge molecular cloud complex is 



located in front of IC 443 (Cornettet al 



1977 



Troia et al.ll2006|). Evidence for a si milar situation 



has also been observed in the ringlike morphology of 3C 391 (IBroean et al 
known to be interacting with a molecular cloud. 



2005), another SNR 




V-tv 

(b) 



18 15 00 17 45 30 15 00 16 45 30 

Fig. 9. (a) A color representation of the n ear-infrared emi ssion observed with 2MASS in the / (in 
blue), H (in green), and K s (in red) bands (R ho et al IboOll The overlaid contours trace the flattest- 



spectrum radio structures of the SNR IC443 between 74 and 330 MHz at a. 



330 



-0.05 and -0.25. 



(b) 330 MHz image of IC 443 showing the locations of the brightest regions of the remnant. 



Towards the interior of the bright eastern shell, the spectral index gradually steepens with po- 
sition in coincidence with a decrease in the intensity of the radio emission. Widely distributed K$ 
band emission is observed in the 2MASS image in this part of the remnant, which was proposed to 
delineate H2 shocked gas from the region interacting with the adjacent molecular cloud. 

In contrast to the excellent agreement between the ionic emitting gas and the flattest spectrum 
features found in the eastern bright limb, no much obvious correspondence is observed in the 
southern part of IC 443 with the exception of a spectral component located at the northern extreme 
of the ridge, near R.A.= 06 h 16 m 45 s , dec.= 22° 35'. This poor IR/radio-spectrum correspondence 
is consistent with the hypothesis proposed before, in which the flat spectrum features in this part of 
IC 443 have a non-thermal origin. 

In the remainder of this section we attempt to infer, using our new measurements at low radio 
frequencies, the physical properties of the area of thermal absorption seen in Fig. [9^ spatially 
coincident with the ionized fine-structure line emitting atoms, i.e. the eastern half of the SNR. At 
radio wavelengths, the emission measure (EM) is given by, 
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EM = 6.086 x 1(T 6 a(T e , v) _i v 2J t v T^' 35 cm" 6 pc, 



(2) 



where a(T e ,v) is the Gaunt factor assumed to be 1, a correct value for the range of astrophysical 
quantities involved in our calculations; v and r v are the frequency measured in MHz and the free- 
free optical depth, respectively, and T e is the electron temperature in K of the intervening ionized 
gas. By measuring the rela tive strengths of the [Fell] lines observed in the near and mid infrared 



emission, 



Rho et al 



(2001) conclude that the emitting region behind the J-type shocks as observed 
in the / and H bands of the 2MASS has a temperature of 12000 K. Although not quantified, the 
authors recognize a large uncertainty associated with this magnitude as a consequence of different 
beam size and possibly different filling factors in their measurements. If we assume that T e is in 
a reasonable range between 8000-12000 K (which includes the temperature as estimated from the 
IR observations) and use the optical depth derived from our radio study particularly for this region 
where the thermal absorption is stronger (T74 ~ 0.3) we obtain an EM between approximately 
2.8 x 10 3 and 5.0 x 10 3 cm -6 pc for the eastern rim. By combining this emission measure with 
the postshock electron density, we can roughly calculate the thickness of the molecular gas layer 
that has been dissociated and io nized by the SNR s hock front . If we assume an ele ctron density of 



500 cm as estimated by 



Fesen & Kirshnerl (119801) and 



Reach & Rhol (|2000) on the basis of 



forbidden [Fell] lines, we conclude that the dissociation and ionization processes took place in a 
thin screen of about 3.4 to 6.0 xlO 16 cm (~0. 01-0. 02 pc). This is a small path compared with the 
transverse dimensions over which thermal absorption is observed, but is a lower limit if the ionized 
gas is clumped. 



6. Comparison with the molecular distribution 

In Fig. [TOl we present a comparis on between our VLA 330 MHz image and new 12 CO (7=1-0) data 
presented by IZhang et al.l (120101) on the basis of new observations carried out with the telescope 
of the Purple Mountain Observatory in China (main beam size of 50" x 54", velocity resolution 
of 0.37 km s , and rms noise level of 0.1-0.3 K at a velocity resolution of ~0.2 km s ). The 
contours superposed on the radio emission depict the CO emission integrated in the range between 
-10 and -1 km s _1 , which includes the systemic velocity of IC 443. As described before, the 
molecular material is preferentially located in the center of the remnant extending in the southeast- 
northwest direction. The spatial distribution of the molecular gas across IC 443 is clearly non- 
uniform. Earlier observations have identified the presence of various clumps of molecular gas with 
broad line widths, as expected from the interaction with t he supernova shock front (clumps labeled 
from A to H in the nomenclature of Dickman et all 1992 ). 



On the basis of the new image at 330 MHz it is possible to recognize details previously un- 
noticed in the spatial comparison between the radio emission and the molecular gas. In Fig. [TO] 
we display the regions where good correlation between radio features and molecular gas distribu- 
tion is observed. Particularly noticeable is Fig. [TOb where it is apparent that the indentation of the 
eastern border in radio occurs near a region where a significant enhancement in the CO emission 
is detected. The molecular complex is transverse to the radio indentation with the maximum of 
the CO emission shifted to the southwest in at least 3' from the border of the SNR (in the region 
of the molecular clump E). It is possible that the singular indentation has formed as the result of 
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the supernova shock front wrapping around a dense clump. Also, Fig.fTOfc shows the presence of a 
concentration in the CO emission around R. A. = 06 h 17 m 16 s , dec.= +22° 25' 40" (molecular clump 
B), that matches a local maximum in the radio emission. The morphological matching between the 
radio synchrotron emission and the molecular gas is especially remarkable at the northern extreme 
of the southern ridge of IC 443 nearR.A.= 06 h 16 m 45 s , dec.= +22° 34' 00" (as shown in Fig. [lOjl), 
where the molecular contours, delineating high density gas in the region of the clump G, are ob- 
served enclosing the bright radio emission. 




OS 19 001 S 30 00 17 30 00 16 30 00 15 30 00 
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Fig. 10. A comparison of the radio continuum emission of IC 443 and the 12 CO (7=1-0) distribution 
in the SNR region. The grayscale representation corresponds to the new VLA 330 MHz image, 
while the superpos ed contours trace CO integrated emission from -10 to -1 km s _1 as taken from 
Zhang et al.l (120101) . Close-ups of three interesting areas are displayed around the center image, 
with the CO contours o verlapping. The white letter in each panel corresponds to the designation of 
the molecular clump bv Dickman et al.1 (119921) . 



We also searched for spectral evidence of shock/cloud interaction. Figure QT]displays an over- 
lay of the 12 CO (7=1-0) integrated emission contours with the radio spectral index distribution 
calculated between 74 and 330 MHz. The CO molecular gas seen in projection onto the plane of 
the sky, overlaps the only flat spectral region observed in the interior of IC 443. Various small com- 
ponents with a spectrum appreciably flatter than the surrounding synchrotron plasma are observed 
distributed nearby or in coincidence with local higher density gas as traced by the 12 CO contours 
along the bright southern ridge and towards the northwest, suggesting that t hese features must be 
region s where strong shocks encountered denser material (as discussed by 



Anderson & Rudnick 



1993b . In the eastern periphery the situation is different. The very flat spectrum component ex- 



tends over an area considerably larger than the molecular cloud, precisely because in this region, as 
shown in Sect. 5, most of the molecules were dissociated and ionized, absorbing the radio emission 
at low frequencies. 
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Fig. 11. A comparison of the spectral index distribution with the molecular emission. The color 
representation corresponds to the spatially resolved spectral index map between 74 and 330 MHz, 
while the overlaid c ontours trace the 12 C0 (7=1-0) emission integrated between -10 and -1 km s -1 
(IZhang et al.ll20101) . The position s where OH(1720 MHz) maser em ission were detected are indi- 



cated with open white diamonds (IHoffman et al.ll2003 



Hewitlll2009l) 



In Fig.QT|we have i dentified with open white diamonds the r egions with OH (172 MHz) maser 



Dubner et al 



d2000J) demonstrated 



emission as detected bv lHewittl d2009l) . In the case of SNR W28,|l 
a clear correspondence between regions of flat spectral index and OH maser emission. We search 
for a similar correlation in IC 443. From Fig. QX]we conclude that there is not a simple association 
between maser locations and spectral index features but rather both flat and steep components are 
observed near masers areas. This fact can be compared with previous results, which suggest that 
maser regions in IC 443 arise from regions with different shock geometry: a shock mostly propa- 
gating towards the line of sight i n the southern OH maser emission and a transverse sho ck in the 

: lHewittll2009l) .Wennd 



westernmost OH emitting region (IClaussen et al 



1997 



Hoffman et al. 2003 



that the OH maser area with transverse shock (R.A.= 06 h 16 m 43 s .6, dec.= +22° 32' 36".7) corre- 
lates well with flat spectral index (ff^ = -0.4), and one of the OH regions where the shock is tan- 
gential to the line of sight (R.A.= 06 h 17 m 29 s .29, dec.= +22° 22' 42" .5) is associated with a steep 
spectrum (ff^ = -0.8). The fact that flat spectrum emission (aif = -0.4) is observed towards 
the easternmost OH maser area with tangential shock (R.A.= 06 h 18 m 3 S .67, dec.= +22° 25' 53".4) 
can be explained because in this region thermal electrons mask the intrinsic SNR spectral effects. 



7. Comparison with TeV emission 

Figure Q~2 ] compares the Te V gamma-ray significance contours obtained from VERITAS obser- 
vations ( lAcciari et al.l 120091) along with the 74/330 MHz spectral index map presented in Fig. [6] 
(Fig. [12li), with the total intensity features of IC 443 at 330 MHz (Fig. [T2J3), and with the 
12 CO (7=1-0) molecular gas distribution (Fig. [T2b). In Fig. [TZb the white plus sign marks the 
position of the source CXOU J061705. 3+222127, while the white cross indicates the centroid of 
the TeV source VER J0616.9+2230. 

Very high energy (VHE) gamma-ray emission roughly extends along the southeast-northwest 
axis, overlapping most of the central region of the remnant. On the basis of the available statistics, 
no correlation can be demonstrated between the radio spectrum or radio features with the TeV 
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emission region. On the contrary, it is remarkable the morphological coincidence between the TeV 
emission and the molecular gas distribution (Fig.fTZt). at least up to the level that the available TeV 
statistic permits to confirm. The interpretation of this striking correspondence is beyond the scope 
of this paper. 
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Fig. 12. Very high energy gamma-rays contours as taken from VERITAS (lAcciari et al. 
perposed with (a) the 74/330 MHz spectral index map presented in Fig. [6] (b) radio continuum 
emission at 330 MHz, and (c) 12 C0 (7-1-0 ) integrated emissi on from -10 to -1 km s~' tracing 
the molecular cloud interaction with IC 443 ( Zhang et alJboiO ). The white plus sign and the white 
cross in (a) mark the positions of the source CXOU J061705. 3+222127 and the centroid of the 
TeV source VER J0616.9+2230, respectively. 



8. Conclusions 

In this work we report on new full-synthesis imaging of the Galactic SNR IC 443 generated from 
multiple-configuration VLA observations at 74 and 330 MHz. These high fidelity images constitute 
the best angular resolution, low frequency radio study published to date on this classic remnant. 
Based on these new data we measured integrated flux densities for this object of S74 mhz = 470 ± 
51 Jy and S330 mhz = 248 + 15 Jy. On the basis of these new total flux density measurements at 74 
and 330 MHz together with previously published values, we recalculated the integrated spectrum 
of IC 443. We have made this analysis, for the first time, taking into account that the spectrum of 
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IC 443 has a reported turnover at the lowest radio frequencies. The fit over a wide frequency range 
produces a radio spectral index Q'!omh^ IHz = -0.39 + 0.01 with free-free thermal optical depth 
T330 ~ 7 x 10" 4 (or T74 ~ 0.02). If the measurement at 10 MHz is excluded, we note that a is 
relatively unchanged from that inferred from a single power law fit (ffJoj^f 112 = -036 ± 0.02). 

Based on the combination of the new images at 74 and 330 MHz, we investigated spectral 
changes with position across IC 443 and related them to the spatial characteristics of the radio 
continuum emission and of the surrounding interstellar medium. For the first time we spatially 

resolved the flattest spectrum region over the SNR (with very flat indices down to al^P 0.05) 

along the brightest rim on the eastern side of the remnant. The radio spectrum steepens towards the 
more diffuse component in the interior of IC 443 with or,? between -0.6 and -0.85, consistent with 
those expected from linear diffusive shock acceleration processes. Finally, another region with non- 
uniform flattening in the spectrum is identified near the southern ridge of IC 443. We investigated 
the connection between these spatial spectral variations and IR and molecular emission distribution. 

The comparison with 2MASS NIR data underscored an impressive coincidence between the 
eastern radio flattest spectrum region and NIR ionic lines which account for most of the infrared 
emission observed towards this part of IC 443. Based on the presence of this IR emission, which 
confirms the existence of a strong J-type shock dissociating molecules and ionizing atoms, we con- 
clude that the most likely explanation for the flattest spectrum observed here, is that it is produced 
by free-free absorption at 74 MHz along the line of sight (with T74 up to ~ 0.3). From the physical 
parameters of the IR emitting gas together with our new radio spectral study, we conclude that the 
observed thermal absorption takes place in a very thin layer (width ~ 0.02 pc) that extends all along 
the eastern border. Such a thin layer would be the product of the dissociating/ionizing action of the 
SNR shock over the adjacent molecular gas. It is important to remark that the low optical depths 
derived for the whole SNR (t 330 ~ 7 x 10~ 4 and r 74 ~ 0.02) indicate that the 74 and 330 MHz 
integrated emission is not significantly attenuated, highlighting that the thermal absorption inferred 
from the spatially resolved spectral index map is a r elatively subtle and localized effect. Our result 
represents only the second case, following 3C 391 ( Brogan et all2005 ). of spatially resolved ther- 
mal absorption delineating the interaction of a SNR/molecular cloud shock boundary. Moreover 
it confirms such phenomena as common and a rich area of investigation for future, low frequency 
studies of Galactic complexes. 

On the other hand, from the molecular studies it is known that the southern ridge, where the sec- 
ond region with flat spectrum was identified, is the site in which the most complex interaction be- 
tween the SNR shock and the external molecular cloud is occurring. We used Zhang et al. (2010)'s 
observations of the 12 CO (7=1-0) molecular gas towards IC 443 to analyze correspondences be- 
tween these data and our new low radio frequency observations. The correlation of the emission 
at 330 MHz with the new 12 CO observations revealed spatial irregularities in the denser molecular 
gas that are well associated with features of bright emission in radio and good agreement between 
the flat radio spectrum region towards the south and molecular emission. We conclude that here 
the flat radio spectrum is predominantly a signature of shock acceleration in a region with strong 
postshock densities and enhanced magnetic fields produced after the interaction of the blast wave 
with dense ambient medium. 

Furthermore, from the comparison of the molecular environment mapped by 12 CO data with 
the local variations of the spectral index, we find evidence that the shocked molecular gas, as 
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illuminated by the OH ( 1 720 MHz) maser emission, is coupled with a flattening in the radio spectral 
index in the locations where the shock is transverse to the line of sight, while a steep spectrum is 
observed in a OH maser region in which the shocks are propagating along the line of sight. 

From the new images, we also analyzed the pulsar wind nebula around the source 
CXOU J061705.3+222127 in IC 443. Based on the new VLA image at 330 MHz we have derived a 
flux density S ^v, H = 0.23+0.05. We have a lso measured the flux density of the PWN in the radio 



330 MHz 

continuum image at 1420 MHz presented bv lLee etalJ d2008h obtaining S^q MH = 0.20 + 0.04. 



From the combination of these data and measurements taken from the literature we derived an 
integrated radio spectrum ar^ ~ 0.0 for the PWN. 

In addition, on the basis of available VERITAS statistics we compared the TeV emission with the 
12 CO distribution, finding an excellent morphological correlation between the high energy emission 
and the distribution of the molecular gas. No correspondence was found, however, between the 
gamma-ray emission as observed by VERITAS and radio spectral or intensity features in radio. A 
future paper will address these findings in connection with the origin of the gamma-ray emission 
detected in IC443. 
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